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ABSTRACT 

We present observations of a large-scale disc of neutral hydrogen (Hi) in the nearby 
Fanaroff & Riley type-I radio galaxy NGC 3801 with the Westerbork Synthesis Radio 
Telescope. The Hi disc (34 kpc in diameter and with Mhi = 1.3 x IO^Mq) is aligned 
with the radio jet axis. This makes NGC 3801 an ideal system for investigating the 
evolution of a small radio source through its host galaxy's cold ISM. The large-scale 
H I disc is perpendicular to a known inner CO disc and dust-lane. We argue that the 
formation history of the large-scale Hi disc is in agreement with earlier speculation 
that NGC 3801 was involved in a past gas-rich galaxy-galaxy merger (although other 
formation histories are discussed). The fact that NGC 3801 is located in an envi- 
ronment of several H I-rich companions, and shows indications of ongoing interaction 
with the nearby companion NGC 3802, strengthens this possibility. The large amounts 
of am bient cold ISM, combined with X-ray results by ICroston. Kraft. &: Hardcastld 
(|2007t) on the presence of over-pressured radio jets and evidence for an obscuring torus, 
are properties that are generally not, or no longer, associated with more evolved FR-I 
radio sources. We do show, however, that the Hi properties of NGC 3801 are com- 
parable to those of a significant fraction of nearby low-power compact radio sources, 
suggesting that studies of NGC 3801 may reveal important insight into a more general 
phase in the evolution of at least a significant fraction of nearby radio galaxies. 

Key words: galaxies: active - galaxies: evolution - galaxies: individual: NGC3801 - 
galaxies: interactions - galaxies: ISM - galaxies: jets 



1 INTRODUCTION 

Over the past decade it has become clear that radio sources 
emanating from super-massive black-holes in the centres of 
galaxies are - rather than merely a by-product of galaxy 
evolution - key players in the formation and evolution 
of galaxies throughout the Universe. It is now fairly well 
established that the properties of Active Galactic Nuclei 
(AGN) are intimately linked with the properties of their 
host galaxies. While empirical evidence has established 
the co-evolution of massive bl ack holes and their host 
galaxies at various rcdshifts ( e.g. Ferra rcsc fc Mcrritt"2000l. 
[Alexander et al. 2005. see also lSilk fc Rees 1998 ). rmmcrical 
simulations demonstrate that AGN-induced outflows are 
important feedback processes for clearing the circum- 
nuclear regions and halting the growth of the super-massive 
black-holes, regulating the correlations found between the 
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black-hole mass and the host galaxy's bulge properties and 
preventing the form ation of too many massive galaxies in the 
early Univer se (e.g. IPi Matte o. Springel. fc H ernauist"200 



Springel, Di Mattco. & Hcrn quist 
20051 . IWagner fc Bicknelil201ll ). 



20051 . .Hopkins ct al 



I 



An important phase of AGN feedback is when syn- 
chrotron jets emanating from the central black-hole region 
start propagating into the host galaxy's inter-stellar medium 
(ISM). At high redshifts, radio sources are observed to be 
aligned with optical continuum f r om w ar m gas and stars 
([Chambers . Milev. fc van Breugej Il987l . iMcCarthv et al.l 
[l987) and the gas in th ese regions is characterised by highly 
disturbed kinematics |McCarthv. Baum. fc SpinradI Il996l. 



Villar-Marti'n. Binette. fc Fosburvi 19991 . Humphrey et alj 



2006f ). This indicates that radio sources vigorously inter 



act with the host galaxy's ISM, entra in and expe l gas and 
induce or quench star formation (e.g. iBicknell et al.,,2000i V 
At intermediate and low redshifts, powerful radio jets are 
know to induce fast ( > 1000km s~^) outflows of i onised 
and neutral gas (e.g. iTadhunteJ Il99ll : IClark et al.l [l998 : 
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Oosterloo et al.' 2000; Morganti et alj 2003 


. l2005allbl. I2OO7I: 


Emonts et alJ 2005i: ,Holt et al. 


2006. 200S 


, I2OIII) and can 


shock-heat the molecular eas 


Guillard et all 201 ll"). exert- 


ing significant feedback on the host galaxy's ISM. 



to be common amo ng more powerful, high-excitation radio- 
loud AGN (see alsolHeckman et"al]|l986l.[Baum et al.ll 19921 . 



While propagating radio jets thus have a profound im- 
pact on their surro unding ISM, the reve rs e is also true (see 
simu l ations by e.g. [Saxton et al. Il2002allbl;rj evakumar et al 



2OO5I : iBicknell fc Su therland 2006; Sutherland & Bicknell 



2003). It is well established that external factors, in par- 
ticular the inter-stellar and inter-galactic medium (ISM 
and IGM), shape the morphological and physical prop- 
erties of radio sources. For example, there is strong evi- 
dence that radio sources in general show a degree of in- 
trinsic asymmetry (apart from orientat ion effects), which 
is most prom inent in compact sources (jSaikia et al ] 1 19951 . 
I2001I . bOOSal V Compact radio sources also show a larger 
degree of polarisation asymmetry than extended sources, 
suggesting a link between the propagatin g radio jets and 
surrounding ISM jSaikia fc Guptal [ioO SbV Also, the im- 
pact of supersonic jets onto the ISM/IGM can create co- 
coons of (shocked) gas and synchrotron emission with 'hot- 
spot' morphologies, while entrainment of gas by slower jets 
can lead to bright ra dio sources with edge- d arkened lobes 
(see e.g. reviews bv [Bridle fc PerlevI 1 19841 : iFerraril 1 19981 : 
IBicknell fc Sutherland 20061 ). These external factors may 



thus contribute to the classical dichotom y in radio sources 
betw een Fanaroff fc Riley type-I (FR-I; iFanaroff fc RilevI 
11974 ) objects, with low-power, sub-relativistic jets that have 
an edge-darkened morphology, and type-II (FR-II) sources, 
with po werful, relativistic jets e nding in a bright hot-spot 
(see e.g. iGopal-Krishna fc Wiita| [2000). 

A particularly interesting object in this respect is 
the nearby {z = 0.0115) radio galaxy NGC 3801. 
NGC 3801 is classified as an FR-I radio galaxy, 
based on its total radio power (P i.4GHz = 2.9 x 



on 

10^^ WHz-i 



— Condon fc Broderickl 1988|) , radio source mor- 
phology (jjenkind 1 19821 : ICroston et al.l l2007l ) and weak 
emission-line spectrum l|Heckman et al.l Il986l). The radio 
source of NGC 3801 is small (~ 11 kpc: |jenkinj|l982h and 
contained well within the optical boundaries of a peculiar 
early-type host galaxy, with a patchy dust feature along 
the major stellar axis and a perpendicular prominent inne r 
dust-lane ijHeckman et al ] | 19861 : IVerdoes Kleiin et al.l[l999l ) , 
as well as a distorted and box-shaped faint outer envelope 
l|Heckman et al.ll 19861 ). The inner dust-lane coincides with a 
kpc-scale CO di sc, orient ed perpendicular to the ra, dio axis 
(|Das et al.ll2005l . see also lOcana Flaauer et aLllioiol ). 

X-ray observations by ICroston et al.l ()2007l ) showed 
evidence for shocked X-ray gas in the central region of 
NGC 3801, created by supersonic, over-pressured jets that 
are not in equilibrium with the hot ISM. These X-ray 
studies also revealed a high absorbing column towards the 
nucleus of NGC 38 01, likely due to an obscuring torus 
ijCroston et al.ll2007l ). In addition, single-dish observations 
sho wed that NGC 3801 is rich in neutral hydrog en (HI) 
gas (|Heckman et al.|[l983l : lDuprie fc Schneideijri996l ). Based 
on its cold gas content and peculiar optical host galaxy 
morphology, it has been suggested that NGC 3801 could 
have been invol v ed in a gas-rich g alaxy-galaxy merger 



jDas et al ] I2OOJ : ICroston et al.l l2007l ). As speculated by 
ICroston et al ] (|2007l ). these are all properties that are gen- 
erally not associated with FR-I sources, but are believed 



iHardcastle. Evans, fc CrostonI 12009 . ICroston et al.l 120081 7 
weU as Sect. 14.4. 1() . 

In this paper we present radio synthesis observations of 
the HI gas in NGC 3801. Our observations reveal that the 
HI gas is distributed in a large-scale disc of cold gas. The 
large-scale H I disc is aligned with the radio axis, indicating 
that the radio jets propagate directly into the cold gas disc. 
Based on our HI results, we discuss the formation history 
of NGC 3801 and investigate further how unique NGC 3801 
is compared to other nearby FR-I radio galaxies. 

Throughout this paper we assume Ho = 71 km 
s"^ Mpc"\ which puts NGC 3801 at a distance of 48.6 Mpc 
and 1 arcsec = 0.24 kpc. 



2 OBSERVATIONS 

HI observations of NGC 3801 were performed with the 
Westerbork Synthesis Radio Telescope (WSRT) on 13 July 
2009 during a single 12h track. A standard calibration was 
applied, using 3C 147 as primary calibrator. Data were re- 
duced and analysed using the MIRIAD and KARMA soft- 
ware. After calibration, the continuum ('channel-0') data 
were separated from the line data by fitting a straight line 
to the line-free channels. A radio continuum map with uni- 
form weighting and beam-size 10.57" x 37.55" (PA -0.4°) 
was produced by Fourier transforming the channel-0 data 
and subsequently cleaning the signal. In the same way, an 
HI cube with robust weig hting +1 l|Briggslll995l ). beam- 
size of 22.47" X 87.52" (PA 0.2°) and channel separation 
of 4.1 km s~^ was created from the line-dataQ We subse- 
quently Banning smoothed the line-data to an effective ve- 
locity resolution of ~ 8 km s~^, resulting in a noise level of 
0.7 mjy bm~^ chan~^. Total intensity images of the H I emis- 
sion were made by using a mask (derived by smoothing the 
positive signal both spatially and kinematically) to extract 
the HI line emission from the original line data-set. The 
data in this paper are corrected for primary beam attenua- 
tion and presented in optical barycentric velocity definition. 



3 RESULTS 

3.1 HI emission in NGC 3801/3802 

Figure [1] show the large-scale HI disc of NGC 3801, which 
has a diameter of 34 kpc. A small fraction of the H I emission- 
line gas in t he dis c was detected with VLA observations by 
iHota et al.l (|2009l ). We derive a total mass of the emission- 
line gas of Mhi = 1.3 X 10® M0, though we note that part of 
the HI disc in our WSRT observations is detected in ab- 
sorption against the radio continuum source (see Fig. [3| 
and not taken into account in this HI mass estimate. The 
HI emission-line gas is aligned with the host galaxy's main 
stellar axis and the east- west dust-lane. The HI disc shows 



^ The large difference in beam-size between north-south and 
east-west direction is a result of the relatively low elevation of 
the source throughout the observations with the east-west array- 
configuration of the WSRT. 
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Figure 1. NGC 3801 and NGC 3802. Left; Total intensity of HI emission (blaclc contours) and 1.4 GHz radio continuum (white contours) 
overlaid into an optical SDSS image of NGC 3801 (see also top-right inset) and NGC 3802. Note that the spatial resolution of the radio 
data in north-south direction is a factor of 4 worse than in east-west direction due to the highly elongated beam (Sect. [2]l. For clarity, 
HI absorption has been omitted from this plot. Contour levels HI: 0.15, 0.21, 0.28, 0.34, 0.45, 0.62, 0.84, 1.1, 1.5, 2.0, 2.5, 3.1, 3.8, 4.5, 
5.2, 6.0, 7.0, 8.0 xlO^Ocm-^, Contours levels radio continuum: 0.04, 0.11, 0.18, 0.25 Jybm"!. Right; Position velocity maps of the HI 
emission (black contours) and HI absorption (grey contours) along the two lines presented in the left image. Contour levels: -9.0, -7.2, 
-5. 6, -4.2, -3.0, -2.0, -1.2 (grey), 1.2, 2.0, 3.0, 4.2 (black) mjybm"^. The cross represents the peak of CO emission from 'clump C found 
bv lDas et all ||2005|) (see Sect.liJl. 



regular rotation with a large velocity spread of Av ~ 600 
km s^^, corresponding to a maximum rotational velocity 
of ~ 300 km s^^ (which is larger th an the rotational ve- 
locity of the optical emission-line gas; iHeckman et al.|[l983l . 
Il985bl ). The HI disc is centred on Vsys = 3452 ± 20 km s"^ 
[z = 0.01151 ± 0.00007), which we argue reflects the sys- 
temic velocity of NGC 3801. Despite the regular rotation, 
Fig. [2] shows that the H I disc is not (yet) fully settled. The 
most prominent asymmetry in the HI distribution occurs 
at the north-western edge of the disc, with gas apparently 
stretching in the direction of the close gas-rich companion 
NGC 3802. 



The close companion NGC 3802 is a disc-galaxy (SO or 
spiral) approximately 33 kpc north of NGC 3801. NGC 3802 
contains Afni ~ 3.9 x lO"^ Mq of HI gas. The HI gas has an 
asymmetric distribution across NGC 3802, peaking in the 
western part of the disc and continuing on the eastern side, 
covering a total velocity range of ~400 km s"'^ (see Fig. 
[T] - right for more details). While the total intensity im- 
age of Fig. [T] (leR) shows a spurious HI bridge in between 
NGC 3801 and the eastern part of the disc of NGC 3802, we 
note that this is merely a projection effect. In fact, Figs. [T] 
(right) and [2] show that there is an apparent bridge or warp 
of H I gas stretching from the north-western part of the H I 
disc of NGC 3801 towards the western part of the disc of 
NGC 3802. This HI feature is too faint to be identifled in 
the total intensity image of Fig. [1] (left ) . The asymmetric H I 
distribution in companion galaxy NGC 3802 and the appar- 
ent HI bridge in between NGC 3801 and NGC 3802 suggests 
that the NGC 3801/3802 system is in ongoing interaction. 



3.2 H I absorption in NGC 3801 

Figure [3] shows the H I absorption detected against the ra- 
dio continuum source in the central part of NGC 3801. 
For the radio continuum we derive a total integrated flux 
of 5*1.4 GHz = 1.03 Jy, corresponding to Pi.4GHz = 2.9 x 
lO^^WHz-^ This is in agreement with earlie r results by 
[Jenkins (198^) and lCondon fc Brodericlj (|l988l ). 

Taking into account that the HI absorption spectra 
shown in Fig.[3]are not entirely mutually independent, there 
are three distinctive features visible. These features are vi- 
sualised in Fig. m and summarised in Table [T] 

The first is a deep and narrow HI absorption ai v — 
3488 ± 20 km s"^ {z = 0.01163 ± 0.00007), which peaks 
at the location of the nucleus. The central velocity of this 
feature is in good agreemen t with the centra l velocity of 
V = 3494 ± 70 km s~^ that W egner et all l|2003l ) derive from 
optical spectroscopy. 

The second is a broad component, covering the central 
region and eastern radio lobe. Figure [1] (middle) suggests 
that this broad absorption component traces H I gas as part 
of the large-scale HI disc, at locations where the gas is in 
front of the radio source|j This also means that the east- 
ern radio jet is most likely the 'far' side of the radio source 
(where the radio continuum is behind the absorbing gas in 
the disc), while the western jet is the 'near' side. 

^ We note that the velocity gradient of the H I absorption as seen 
in Fig. [T] could be slightly steeper than that of the H I emission. 
This may occur because the spatial resolution of our observations 
is coarser than the underlying radio continuum structure, hence 
the H I absorption is less affected by beam-dilution than H I emis- 
sion and could trace the more inner regions of the H I disc. 
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Figure 2. Channel maps of the HI emission in NGC 3801, binned to a velocity resolution of 43 km s ^. The velocity range of the 
HI absorption is largely omitted from this plot. Contour levels: -3.2, -2.6, -2.0, -1.6, -1.2, -0.8 (grey), 0.8, 1.2, 1.6, 2.0, 2.6, 3.2 (black) 
mjy bm~^ . 



The third feature is a narrow absorption against the 
western jet, which is blueshifted by ~ 50 km s^^ with re- 
spect to Vsys- Fig. [T] (middle) shows that also this HI com- 
ponent is most likely part of the large-scale, edge-on HI 
disc that embeds the small radio source, although it cannot 
be ruled out that this blueshifted HI absorption represents 
an outflow of neutral gas driven by the propagating radio 
jets. Figure [3] also shows indications of a narrow and much 
weaker HI absorption component against the eastern radio 
jet in NGC 3801, around v~ 3550 km s'V 

Higher resolution H I absorption observations are essen- 
tial for studying the nature of the various absorption fea- 
tures in more det ail. For this, we refer to ongoing work by 
iHota et all l|2009l) . 



3.3 H I environment 

Figure [S] shows the HI environment of NGC 3801, with five 
other galaxies detected in HI emission. Table [2] summarises 
the HI properties of these systems. The bulk of the HI in 



these companion galaxies is distributed in disc-like struc- 
tures. 



We also find two 'clouds' of HI gas that are not associ- 
ated with an optical counterpart in the SDSS image of Fig. 
[S] One of them, cloud A, has a velocity that coincides with 
that of the eastern part of the HI disc of NGC 3801, as well 
as companion 3, and is located 83 kpc south of NGC 3801 
(note that the velocity of cloud A is significantly offset from 
that of the much nearer companion 4; see Fig. [5] - left) . 
Cloud B coincides in both velocity and location with the 
H I-rich companion 6 (NGC 3806) and is most likely gas that 
is either stripped off, or being accreted onto, the HI disc of 
NGC 3806. If these clouds of HI gas are tidal in origin, and 
given that the total H I mass of both clouds (4 — 5 x 10^ M0) 
is similar to that of NGC 3801 's close companion NGC 3802, 
this indicates that past as well as ongoing interactions could 
have a profound influence on the evolution of this group of 
galaxies. 
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Figure 3. HI absorption properties of NGC 3801. The middle plot shows a zoom-in of the optical host galaxy, HI emission and radio 
continuum from Fig. [T] The zoom-windows show the H I profile at various locations along the H I disc. 



4 DISCUSSION 



4.1 Large-scale HI disc 

Our HI results show that the bulk of t he HI gas de- 
tected in earlier sinsle-dish ob servations (jHeckman et al.l 
1 19831 : iDuprie fc Schneideilll996t ) is arranged in an H I disc 
with a diameter of 34 kpc, aligned along the host galaxy's 
major stellar axis and prominent east-west dust-lane. We 
therefore argue that this large-scale HI disc traces the 
major axis of the host galaxy. The peak HI surface den- 
sity in the disc is Ehi ~ 2.1 Mq pc~^. This is below (or 



at best close to) the critical HI surface densi ty for star 



1989*: 'van der Hulst et al.lll993 


: iMartin & KennicuttI l200ll: 


Bigicl ct al. 2008; Leroy et al.l 
models fe.g. ISchavell2004). It is 


20081). as also predicted bv 
therefore unlikely that wide- 



spread star formation is happening in the disc (although star 
formation may occur in regions of higher local density). 

The rotation of the large-scale HI disc allows us to 
estimate the dynamical mass of NGC 3801 to be Mdyn = 
(RvVG) • sin-^i = 3.6 x 10^^ Mq (with R= 17kpc the ra- 
dius of the HI disc, v = 300km s~^ its rotation velocity 
and assuming the HI disc is viewed edge-on, i.e. i — 90°). 
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Figure 4. Gaussian fits to the HI absorption profiles of Fig. |3] The black solid line shows the total fit, while the red lines show the 
individual Gaussian components. The numbers correspond to the absorption features as listed in Table [l] The actual values in Table [T] 
are derived from the Gaussian profiles with solid red lines and thick numbers in this plot. 



Table 1. HI absorption properties 

Feature Vcentro FWHM 
# (km s-i) (km s-i) 



Nhi 
(xl02" cm-2) 



3488 
3564 
3399 



33 
136 
42* 



0.035 
0.028 
0.014 



2.1 
7.0 
1.1* 



Notes - Values have been derived from the Gaussian fits to the 
profiles, as shown in Fig. |4] Vccntrc is the central velocity of the 
absorption feature (error ±20 km s"-"^); FWHM is the full width 
at half the maximum intensity; t is the maximum optical depth 
that the feature displays across the radio continuum and is 
calculated through = 1 — -^ats. (with Scont the underlying 
21cm radio continuum fiux at the limit spatial resolution of our 
observations); A^hi is the derived HI column density (assuming 
Tlpin = 100 K). 

* derived from both Gaussian components that comprise feature 
3 in Fig. |4] 



Giveri that Lb = 2.3 x 1O^°L0 (|Heckman. Cartv. fc BothunI 
Il985al . corrected for Ho = 71km s-^Mpc-^), NGC 3801 
has Mdyn/is ~ 15, which indicates that NGC 3801 con- 
tains a significant dark-matter halo (in agreement with 
M/ L ratios observed in other Hl-rich early-type galajc- 
ies; Bertola et al. 19931. iFranx. van Gorkom. fc de Zeeuwl 



11994 iMorganti et al.lll99llWeiimans et al.ll2008l ). Consid- 
ering only the HI emission, we find Mhi/Lb = 0.06 for 
NGC 38010 which is within the ran ge of values observed 
for o ther Hl-rich early- type g alajcies (lOosterloo et al.ll2007l . 
and polar-ring galaxies (|Bettoni et al.ll200ll ). 
The radio source is aligned roughly along the major axis 
of the H I disc, indicating that the radio jets are propagating 
directly into the dis c. Perpend ic ular to the large-scale HI 
disc and radio axis, iDas et al.l (|2005l ) detected a 4.6 kpc 
wide CO disc, which follow s a prominen t inner du s t-lane 
(|Verdoes Kleiin et al.l ll999"). We refer to iDas et all l|2005l , 
their figure 5) for an excellent illustrative image of the inner 
region of NGC 3801 (which shows contour plots of the CO(l- 
0) and VLA 20cm radio continuum overlaid onto an HST 
image of the inner dust features). 



Note that this estimate does not take into account the H I gas 
detected in absorption, therefore the true value of Mhi/^b is 



somewhat larger. 



iDas et al.l (|2005l ) also detected a CO clump that is not 
part of the inner CO disc and aligned along the EW dust- 
lane, with a projected distance of about 2 kpc from the nu- 
cleus. They dismiss the likelihood that this CO clump is part 
of a second disc, given that it would require the geometry of 
the inner and outer disc to be extremely rare (roughly one 
in a thousand), namely such that both discs would have to 
be edge-on, within 10° of the line-of-sight and perpendicular 
to each other to within 10°. Our HI results, however, show 
that NGC 3801 does in fact satisfy the conditions of two 
perpendicular, edge-on discs (an inner kpc-scale CO disc 
and a large-scale outer HI disc). Figure [T] shows that the 
CO clump detected by iDas et al] (2005) likely represents 
cold molecular gas that is part of the large-scale outer disc. 
Given that two discs with intersecting orbits are not stable, 
it is likely that the large-scale H I disc has a central gap and 
thus may not continue all the way down to the kpc-scale 
central region. 

Our HI results (Fig. [1]- middle and Fig. [2]), however, 
show that there is neutral gas in the central region, which 
has a steep velocity gradient relative to the bulk of the gas 
in the large-scale H I disc. This H I component is most clear 
in Fig. [2] at v = 3165 km s"^ (panel 1), v = 3251 km s"^ 
(panel 3) and v = 3728 km s~^ (panel 11), of which the lat- 
ter two correspond to the velocity of the two clumps of CO 
^s that represent the central molecular disk in D as et al.l 
(2005). This could indicate an HI counterpart to the central 
CO disc, although the spatial resolution of our H I observa- 
tions is too poor to verify this. 

From an HI perspective, NGC 3801 is not unique. Simi- 
lar to NGC 3801, examples of other early-type systems with 
a perpendicular inner and outer disc are known to exist, 
possibly related to the tri-axial nature of the mass dist ribu- 
tion. These include NGC 5266 jMorganti et al.lll997l ) and 
NGC 2685 (where the two discs/rings are likely a projec- 
tion effect of an extremely warped disc: iJozsa et al.i ,2009). 
NGC 3801 also bea rs resemblance to pol ar-ring galaxies, 
such as NGC4650A l Arnaboldi et al.| ll997l ) and UGC 9796 
l|Cox. Sparke. fc van Moorscl 2006. ) . 



4.2 Formation history and AGN triggering 

Based on its peculiar optical morphology and CO content 
in the inner few kpc, it has been suggested that NGC 3801 
could have been involved in a gas-rich galaxy-galaxy merger 
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Right Ascension {J2000) 



Figure 5. HI environment of NGC 3801. Left; Contours of the total intensity HI emission (black) overlaid onto an optical SDSS image. 
Contour levels HI: 0.15, 0.21, 0.28, 0.34, 0.45, 0.62, 0.84, 1.1, 1.5, 2.0, 2.5, 3.1, 3.8, 4.5, 5.2, 6.0, 7.0, 8.0 xlO^Ocm-^. The 1.4GHz radio 
continuum source is represented by white contours in the very central region of NGC 3801 (levels: 0.04, 0.11, 0.18, 0.25 Jybm~^). For 
clarity, HI absorption has been omitted from this plot. Right: Position- velocity plots of the same HI data cube, represented as total 
intensity of the HI emission integrated along the RA and dec axes. 



Table 2. HI companions 



# 


Name 


D (kpc) 


Vsys (km s 1) 


Av 


Mm (M0 


1 


NGC 3801 




3552 ± 20 


607 


1.3 X 10^ 


2 


NGC 3802 


33 


3250 ± 40 


398 


3.9 X lO'^ 


3 


J114011. 27+174846.6 


75 


3674 ± 20 


103 


3.1 X 10* 


4 


KUG 1137+179 


49 


3858 ± 40 


308 


2.0 X 10^ 


5 


NGC 3790 


102 


3373 ± 40 


380 


6.3 X 10* 


6 


NGC 3806 


118 


3486 ± 20 


180 


4.0 X 10^ 




'Cloud A' 


83 


3657 ± 20 




3.8 X lO'^ 




'Cloud B' 


132 


3496 ± 20 




5.0 X lO'^ 



Notes - D is the projected distance to NGC 3801; Vsys is the central optical barycentric velocity derived from our HI data; Av is the 
velocity range covered by the HI emission; Mhi is the total HI mass detected in emission. 



(e.R. iHeckman et all Il986l : iDas et al.l l2005l : ICroston et al] 
|2007| ). Our HI results are in agreement with such a sce- 
nario, given that the formation of a large-scale H I disc can 
be the natural outcome of a galaxy merge r over significan t 
time-scales (e.g. lEmonts et al.„ 2006 2008: Seira et al.ll200d) . 
Simulations of merging galaxies bv lDi Matteo et al.l ( 20071 ) 
show that gaseous and stellar discs in the progenitor gala^xies 
can be tidally disrupted, depending on the orbital parame- 
ters of the merging galaxies (in partic ular dir e ct enc ounters 
can create large gaseous tidal-tails). iBarned l|2002l ) shows 
that part of this expelled material can be re-accreted back 
onto the newly formed host galaxy and settle into a large- 
scale gas disc on time-scales of one to several Gyr. The bulk 
of the HI gas in the disc of NGC 3801 shows regular rota- 
tion, which means that the gas must have had enough time 
to settle into the observed disc-like structure. As a rough es- 



timate of the merger time-scale, we assume that this process 
must have lasted at least one orbital period, or ^ 3.5 x 10* 

yr- 

The age of the r adio source is signif icantly less than 
this (s; 2.4 X 10'^ vr: ICroston et alJl2007l ). While there is 
no evidence for a direct causal connection between the pos- 
sible merger event and the triggering of the radio source, 
there would have been a significant time-delay between the 
start of the merger and the onset of the current episode 
of the radio-AGN activity. Similar time-delays between 
merger (and associated starburst episode) and radio-AGN 
activity have been observed i n other ne arby radio galaxies 
(e.g.lEvans ct al. 1999; Tad hunter et al.| [2005: Emonts ct al] 
I2006n2008i : iLabiano et al.ll2008l ) 

We note, however, that our HI results alone do not 
provide conclusive evidence that the large-scale HI disc in 
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NGC 3801 was formed as a result of a major merger. Minor 
mergers or interactions with satellite galaxies (similar to the 
ongoing interaction between NGC 3801 and NGC 3802) may 
also have played a role in both shaping the H I disc as well 
as triggering nuclear activity in NGC 3801. Given that the 
bulk of the HI gas that we find in our data cube is associ- 
ated with companion galaxies that are all within v ~ 300 
km s~^ from the central velocity of NGC 3801, it is possible 
that the HI gas in the disc of NGC 3801 was delivered by 
former gas-rich companions ove r the course of seve ral Gyr 
(similar to what is discussed bv lStruve et al.ll2010al for the 
large HI disc in radio galaxy NGC 1167). 

Another possibility is that accretion of cold gas played 
a major role in t he fo rmation of the HI disc. Simulations 
bv lMaccio et al.l (|2006l ) show that extended rings of gas in 
polar-ring galaxies may form naturally through the accretion 
of cold gas falling in along large filamentary structures, while 
[Bournaud & Combos (2003) argue that it is more likely for 
polar-ring galaxies to acquire their gaseous discs through 
tidal accretion from a gas-rich donor galaxy rather than a 
violent merger. A similar mechanism may have formed the 
large-scale HI disc in NGC 3801. 



4.3 Radio source evolution 

As mentioned above, the alignment of the small radio source 
in NGC 3801 with the large-scale outer HI disc makes 
NGC 3801 an ideal nearby radio galaxy for studying the 
early stage of radio-source evolution and its interaction with 
the surrounding co l d ISM . 

ICroston et al.l (|2007l ) find evidence for shocked heated 
shells of hot gas associated with the radio source in 
NGC 3801, caused by supersonic, over-pressured radio jets. 
Their findings indicate that the radio source in NGC 3801 
is heavily interacting with the surrounding ISM. In pro- 
ceedings b a sed o n high-resolution VLA HI observations, 
iHota et al.1 (|2009l ) report the presence of a broad, faint, 
blueshifted absorption and an H I clump associated with the 
shocked shell around the eastern lobe, possibly reflecting 
a gas outflow. Our lower resolution HI data do not show 
clear evidence for a broad, blueshifted absorption, but such 
a faint H I absorption feature could be flUed in with H I emis- 
sion within the relatively large beam of our observations. We 
do detect an apparent blueshifted narrow absorption com- 
ponent against the western radio lobe, but as discussed in 
Sect. 13.21 this most likely represents HI gas as part of the 
large-scale HI disc. 

Similar to NGC 3801, there are other known radio 
galaxies in which the radio source propagates into an 
outer H I disc. This includes Coma A, where the radio 
lobes ionise part of the disc (|Morganti et al.l l2002t ) and 
IC 5063, where the rad i o jets produce fa,st out flows of H I gas 
l|Morganti et al.1 1 19981 : lOosterloo et al.1 12OO0I '). Particularly 
interesting is also the similarity with the disc-dominated 
radio galaxy B2 0722-1-30, which hosts a small FR-I radio 
source (similar to NGC 3801) that propagates close to the 
plane of the disc, in the direction of a heavily interacting 
Hl-rich galaxy pair l|Emonts et al.ll2009l ). 



4.4 H I-rich radio galaxies: the general picture 

In general, roughly 10% of early-type galaxies outside clus- 
ters contain large amounts of HI gas (Mhi ^ 10^ Mp ), often 
distributed in large-scale d i sc- or ring-like structu res (ISadlej 
2001'; 'Morganti ct al.l l2006l : lOosterloo et"al]|2007l . l20loM 1ln 
Emonts ct al. (2010), we showed that the overall HI proper- 
ties of nearby, non-cluster low-power radio galaxies are sim- 
ilar to those of radio-quiet early-type galaxies, with large- 
scale H I discs detected in a significant fraction of radio- loud 
systems. In addition, we showed that there is a clear segrega- 
tion in HI content with radio size, in the sense that massive, 
large-scale HI discs are only found in the host galaxies of 
compact radio sources, while none of the more extended FR- 
I radio sources stud i ed co ntain similar amounts of HI gas 
(|Emonts et al.ll2007l . I2OIOI ) . The ra dio source in NGC 3801 
is ~ 11 kpc in diameter (| Jenkins! J 9821. As illustrated in 
Fig. |6] the HI properties of NGC 3801 are thus similar to 
those of a significant fraction of nearby low-power compact 
radio sources. 

It is likely that in these systems the early evolution of 
the radio source is intimately linked with the surrounding 
cold ISM, and that either confinement by the dense ISM 
or inefficient fuelling (through dumpiness of the cold ISM 
and associated non-steady accretion) keeps the radio source 
compact (see discussions in lEmonts et al.ll2007 ^ 2010l). Anal- 
ternative scenario, recentlv raised bv lMorganti et al.l (|201lD 
in a study of the powerful compact radio source PKS 1814- 
637, suggests that the radio emission in these H I-rich com- 
pact radio sources may be boosted, at least temporarily, 
by the interaction of the radio jets with the rich ISM, im- 
plying that intrinsically these compact source s may refiect a 
link t o radio- AGN in Seyfert galaxies (see also lEmonts et all 
I2OIOI , where we speculate that the transition from radio- 
AGN in Seyferts to the brighter compact radio sources to 
classical FR-Is may perhaps be linked to a transition from 
well-defined star-forming H I discs to low surface brightness 
H I discs to H I-poor host galaxies) . 



4.4.1 NGC 3801: a classical FR-I? 

NGC 3801 is classified as an FR-I radio galaxy, based on 
its total radio power (- Pi.4GHz = 2.9 x 10^^ W Hz~ ^), ra- 
dio source morphology (|jenkinslll982l: ICroston et al. 120071 

1986 1. 



and weak emission-line spectr um (iHeckman et al 



However, X-ray observations bv lCroston et al. I (l2007l l show 



that NGC 3801 has supersonic, over -pressured radio jets , 
not commonly found in FR-I sources (jCroston et al.l l2008ll. 
These X-ray studies also reveal a high absorbing column 
towards the nucleus (possibly representing an obscuring 
torus), somet hing that is generally not associated wit h 
FR-I sources (iHardcastle et al.ll2006l : ICroston et al] l2007l l . 
ICroston et al.l ( 2007 1 suggest that these properties may 
more closely resemble those of high-power, high-excitation 
FR-II sources, which are believed to be post-merger systems 
that are fuelled by cold gas driven into the central region to 
form a radiatively-efficient classical accretion disc (contrary 
to low-power/low-excitation FR-Is, which are believed to be 



■* Idi Serego Aliehieri et"all ll2007l l and lOosterloo et al.l 1I2OIOI I 

show that the H I detection rate of early- type galaxies in the Virgo 
Cluster is dramatically lower. 
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Centaurus A (see also Kraft et alj 20031: van Gorkqm et aD 



10 100 
radio size (Icpc) 



1000 



Figure 6. Adapted from lEmonts et al] 1I2OIC1I '). Shown is the HI 
mass content (or upper limit) plotted against the total linear ex- 
tent of the radio source for a complete sample of nearby non- 
cluster, low-power (FR -I-type) radio galaxie s. T he square rep- 
resent s NGC 3801. See lEmonts et al.l | |2007^ and lEmonts et al.l 
l|201(]| ) for details. For completeness, radio host galaxies des- 
ignated as elliptical (E) in the NASA Extra-galactic Database 
(NED) are plotted in grey, while non-ellipticals (SO and merger 
systems) are black. While this reflects the long known idea that 
the relative HI fractio n seems to rise from ellipticals to SO and 
later-type systems (e.g. lWardle fc Kn app 1986), it also shows that 
the H I-rich compact radio sources in this sample are hosted by 
galax ies with a variety of optical morphologies (see lEmonts et aU 
I2OI0I . for detailed optical imaging). 



fuelled mainly through a quasi-spherical Bondi accretion of 
circum-galactic hot gas directly onto th e nucleus; see e.g; 
Heckman et al.l Il986l iBaum et al 1992, Hard castle et al 



2006 



Hardcastle. Evans, fc Croston. .2007.. .Baldi fc Capetti 



200alEmonts et al.ll20ld . lRamos Almeida et al.ll2011^ ■ This 
scenario is in agreement with our HI results, which show 
that NGC 3801 is likely a post-merger system with large 
amounts of cold gasQ 

However, as we saw in Sect. 14.41 the presence of a 
large-scale HI disc is not unusual for low-power compact 
radio sources like NGC 3801. We speculate that the pos- 
sible link to radio- loud Seyferts (discussed in Sect. 14. 4p 
could perhaps also explain some of the unusual char- 
acteristics of NGC 3801. Many Se yfert-AGN also show 
heavy obscuration in X-rays (e . g. [Maiolino et al.l [l99^ . 
iRisaiiti. Maiolino. fc Salvatil Il999l . ICappi et all 12009 ) and 
can be found in disc-dominated, int eracting galaxies that 
contain large amounts of HI gas (e.g. IKuo et aLll2008l ). 

Regardless of the exact nature of the radio- loud AGN, 
the combined X-ray and HI results suggests that NGC 3801 
is a system with over-pressured radio jets that are embedded 
in (and interacting with) a rich, ambient cold ISM that has 
been deposited after a gas-rich merger; features not or no 
longer pre sent in more e v olved FR-I radio sources. As dis- 
cussed by ICroston et al.l (|2007l ) , this resembles the case of 



^ Our H I results do not show evidence that cold gas is be- 
ing accreted onto the central region (see Section 13.211 , but note 
that the accretion rate necessary to sustain AGN activity could 
be so low that it is w ell beyond our observing capabilities 
llvan Gorkom et al.lllQSgh . 



I1990I ; ISchiminovich et al.lll994l ; IStruve et al.ll2010bl ). X-rav 
observations of other H I-rich compact FR-I radio sources 
may reveal how unique NGC 3801 (and Cen A) are in the 
evolution of classical radio sources. 



5 CONCLUSIONS 

We have shown that the H I gas in NGC 3801 - detected 
in e arlier single-dish observatio ns bv lHeckman et al] ([1983) 
and iDuprie fc SchneideJ l|l996l ) - is distributed in a regu- 
larly rotating large-scale HI disc. We also detect HI gas in 
the close companion NGC 3802 (which is in apparent on- 
going interaction with NGC 3801) as well as several other 
galaxies and gas clouds in the immediate environment of 
NGC 3801. The small radio source in NGC 3801 is propagat- 
ing directly into the large-scale HI disc, making NGC 3801 
an important system for investigating in detail the evolu- 
tion of a small radio source through its host galaxy's cold 
ISM. Because of the similarity between the HI properties 
of NGC 3801 and those of other nearby low-power compact 
radio sources, studying NGC 3801 may reveal important in- 
sight into the early evolution of at least a significant fraction 
of low-power radio sources. 
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